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Abstract 



In the present article we consider a capillary compressible system introduced by 
C. Rolide after works of Bandon, Lin and Rogers, called the order-parameter model, 
and whose aim is to reduce the numerical difBculties that one encounters in the 
case of the classical local Korteweg system (involving derivatives of order three) or 
the non-local system (also introduced by Rohde after works of Van der Waals, and 
which involves a convolution operator). We prove that this system has a unique 
global solution for initial data close to an equilibrium and we precisely study the 
convergence of this solution towards the local Korteweg model. 

1 Introduction 

1.1 Presentation of the systems 

In the mathematical study of liquid-vapour mixture, Gibbs first modelled phase transi- 
tions thanks to the minimization of an energy functional with a nonconvex energy density 
(see [17]). The phases are separated by an hypersurface and there are mainly two ways 
to describe it: either we consider that the interface behaves like a discontinuity for the 
fluid parameters (this is the Sharp Interface model), either we consider that between the 
phases lies a thin region of continuous transition (this is the Diffuse Interface approach, 
where the phase changes are seen through the variations of the density and which is much 
simpler numerically) . Unfortunately the basic models provide an infinite number of solu- 
tions (few of them being physically relevant) and this is why authors tried to penalize the 
high variations of the density (with capillary terms related to surface tension) in order 
to select the physically correct solutions. 

In the present paper, we are interested in the local and non-local Korteweg systems 
(in the diffuse interface model). These systems are based upon the compressible Navier- 
Stokes system with a Van der Waals state law for ideal fluids, and endowed with a 
capillary tensor. 

Let us recall that the local model was introduced by Korteweg and renewed by Dunn 
and Serrin (see [16]) and the non-local model was introduced by Van der Waals and 
renewed by F. Coquel, D. Diehl, C. Merkle and C. Rohde. For an in-depth presentation 
of the capillary models, we refer to [30] and [11]). 
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Let p and u denote the density and the velocity of a compressible viscous fluid {p is 
a non- negative function and u is a vector- valued function defined on W^). We denote by 
A the following diffusion operator 

Au = pAu + (A + p)'Vd[vu, with ^ > and v = X + 2p > 0. 

The Navier-Stokes equations for compressible fluids endowed with internal capillarity 
read: 

J dtp + div (pu) = 0, 

\ dt{pu) + div {pu ®u) - Au + V{P{p)) = KpVD[p]. 

The capillary coefficient k may depend on p but in this article it is chosen constant. In 
the local Korteweg system (NSK), the capillary term D[p] is given by (see [16]): 

D[p] = Ap, 

and, in the non-local Korteweg system {NSRW) (see [29], [11], and [34]), if cf) is an 
interaction potential which satisfies the following conditions 

{\.\ + \.\'^)4>{.) € L\R''), / (l){x)dx = l, even, and > 0, (1.1) 



then D[p] is the non-local term given by: 

D[p] = (j)*p- p. 

Comparing the Fourier transform of the capillary terms, we have {4){Cj — ^)p{C} iii the 
non-local model, and — |Cl^p(C) ™ the local model so that a natural question is to study 
the closedness of the solutions of these models when is formally "close" to 1 — |.^|^. 
For this, we introduced in [7] a specific interaction potential and considered the following 
non-local system: 



{NSRWe 
with 



dtPe + d\w{peUe) = 0, 



dtipeUe) + div {peU ® Ue) - AUe + V{P{pe)) = Pe^V {(t>e * Pe - Ps), 



l0(^) With 0(x) = ^ 



MO 



£2 



For a fixed ^ the Fourier transform of (j)^ is (j)£{S,) = e ^ 151 ^ and when e is small, 

is close to —[Cla- 
using energy methods, we proved in [7] that this system has a unique global strong 
solution for initial data close to an equilibrium state. The functional setting are classical 
and hybrid Besov spaces (taylored to the capillary term). We also obtained that when 
the small parameter e goes to zero, the solution tends to the corresponding solution 
of the local Korteweg system and we obtained a rate of convergence in terms of e. In 
[9], we provided by Lagrangian methods more precise a priori estimates giving a better 
understanding of the convergence and the hybrid Besov setting in terms of the linear 
Fourier structures. 
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Though, these models are not completely satisfying. On one hand, recalling the re- 
sults from [12], [2] (compressible Navier-Stokes system), [14], [8] (local Korteweg model), 
and [19], [7] (non-local Korteweg model), we observe that the density in the local cap- 
illary model is far more regular than in the non-local model where it shares the same 
frequency structure as in the compressible Navier-Stokes model (heat regularization in 
low frequencies and only a damping in high frequencies). 

On the other hand, from a numerical point of view the local model is difficult to 
handle because the capillary term contains third-order derivatives. The non-local model 
also presents difficulties in numerical studies: even if the capillary term only contains 
derivatives of order one, it involves a convolution operator, whose numerical difficulty is 
comparable. 

For this reason, C.Rohde presented in [31] a new model, called the order-parameter 
model, and inspired by the work of D. Brandon, T. Lin and R. C. Rogers in [5]. 

This new system consists in introducing in the capillary term q:^V(c — p) a new 
variable c called the " order parameter" , which is coupled to the density via the following 
relation related to the Euler-Lagrange equation from the variational approach (a controls 
the coupling between p and c): 

e^Ac + a^{p-c) = 
so that the new system he considers is the following: 

{ dt Pa + div (paUa) = 0, 
aUa ^ Ua) ~ ■A.Ua + ^(^(Pa)) — '^'^ Pa'^ {pa Pa)i 
e'^Aca + a'^ipa - Ca) = 0. 

As emphasized by C. Rohde, from a numerical point of view this system is much more 
interesting because now we only have one derivative in the capillary tensor (which is 
local), and the additionnal equation for the order parameter is a simple linear elliptic 
equation that can be easily and numerically fast solved at least when the mesh is fixed. 
Moreover as we will see later, as for the previous non-local capillary model, this system 
has the same frequency structure as the classical Navier-Stokes model. 

In [31], C. Rohde proves (for £ = X = p = lm the two-dimensionnal case) that the 
system has a unique local classical solution: 

Theorem 1 ([31]) Assume that the initial data {pQ,U()) is independant of a > with 
uq G iJ^(]R^), pq > and po — p G H^(M?) for some constant p. Let cq the solution of 
the elliptic problem —Acq + o^cq = cP' po- There exists a constant T=k > such that the 
initial-value problem (NSOPa) has a unique solution defined on [0,T^[ satisfying: 

Pa-p,UaeL'^{0,T,;H\R^)), > 0, 

Ca- pe L'=^{0,T^;H^{M.'^)). 
Moreover, for all t G [0, T^, [ we have 

lim \\pa{t, .) - Ca{t, .)IIl2(R2) = 0. 

In this paper, C. Rohde also conjectures that when the coupling constant a goes to 
infinity, these solutions converge to the solution of the local Korteweg model. 
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1.2 Statement of the results 

In the present article, following what we did in the whole space for the non-local system 
(see [7]) and using lagrangian methods from [9], we will prove that under smallness 
conditions, and with less regular initial data, the system has global strong solutions in 
the following critical spaces (we refer to the appendix for more details on Besov spaces 
and hybrid spaces). We also prove the above conjectured convergence and give an explicit 
rate of convergence with respect to a. 

Definition 1 The space is the set of functions {q, c, u) in 

(a(M+, n Bl,) n Li(R+, s^+i'^-i n K+''')y x (c,(r+, bI;,^) n l1(m+, 

endowed with the norm \\{q, c,u)\\fs = \\{q,c,u)\\psi^^-j where for all t we denote (recall 
that vq = min(/u, v) ) 

de / 

+ z^ohlljiB^+i + u\\q\\ii^s+i,s-i + i^^\\q\\i^^ ^s+2,s + z/||c||2;i^s+i,.-i + z/2||c||2;i^.+2,s (1.2) 

Theorem 2 Let q > and assume min(/i, 2/i + A) > 0. There exist two positive con- 
stants rjop and C only depending on d, fi, X, k, and P'ij)) such that for all rj < rjop, if 

d -j^ ^ d_ . d ^ 

Pq —'p ^ B21 n i?2^i, no G -62^1 ' '^0 Is defined by — Acq + a^co = a^po and 

IIpo - pII . . d + lluoll . d_i < ?7 

then system (NSOPa) has a unique global solution (p^, Cq, Uq,) with {p^ — p,Ca — p,Ua) S 

d 

Fa such that: 

\\iPa-p,Ca-p,Ua)\\ d < C° C {\\po - p\\ d_, d +||no||.d_J. 

nB^^ Bl^ 

Moreover we have the global in time results: 

+T^\\Ca-Pa\\- d > 0, 

L°°(R+,BJ^ ) L°°(IR+,_B2 j) 

II d, +y'^\\Ca-Pa\\ d < C°a"^. 

Bh ) LHR+,Bl,) 




The following result deals with the convergence in a: when the initial data are small 
enough (so that we have global solutions for (NSK) and (NSOPa)) the solution of 
(NSOPa) goes to the solution of (NSK) when a goes to infinity. 

Theorem 3 With the same assumptions as before, there exists < rjQ < mm{riK,VOp) 
such that for all r] < r]Q, if 

IIpo - pII . d 1 d + \\uo\\ . d_i < r], 

BL ^Bl, Bl, 
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then systems (NSK) and (NSOPa) both have global solutions and \\ {pa — p,Ca — p, Ua — 
u)\\ d goes to zero as a goes to infinity. Moreover, with the same notations as before, 

^ a 

there exists a constant C = C{rj,K,'p,P'{l)) > such that for all h g]0, 1[ (if d = 2) or 
h g]0, 1] (ifd> 3) 

\\{Pa - P,Ca - P,Ua - U)\\ d_^<Ca~^, 

Remark 1 We can assume that e = 1 without loss of generality. If not we just have to 
replace k by and a by a/e. 

Remark 2 As the order parameter Ca goes to pa, we formally get that when a goes to 
infinity, the capillary term goes to npVAp. 

1.3 outline of the paper 

The article is structured the following way: section 2 is devoted to the proof of theorem 
2. We first introduce an interaction potential (pa that allows us to rewrite the system 
into a non-local shape. As we want precise estimates we follow the methods from [9]: 
we first obtain estimates on the linearized system and then on the advected linear sys- 
tem thanks to a Lagrangian change of variable. The rest of the proof is classical, we 
define approximated solutions thanks to the Friedrichs' scheme and obtain existence and 
uniqueness like in [7]. In section 3 we prove theorem 3 and in the appendix, we first re- 
call basic properties of Besov spaces, then we provide estimates for the flow of a smooth 
vectorfield. The last part of the appendix is devoted to Bessel functions that are needed 
for the expression of our new interaction potential. 

2 Proof of theorem 2 
2.1 Interaction potential 

As announced in the introduction, we first rewrite the system in a non-local shape. Let 
us focus on the last equation, we can write that (for more clarity we drop the subscripts 
with a): 

-e^A{p - c) + a^(/j - c) = -e'^Ap 

which leads to: 

a'^ip - c) = -a'^i-A + %y'^Ap = -{-^ + ^)"^A/9 
So that in Fourier variable: 

Then up to choose k = and replace a by a/e, from now on we assume that e = 1 and 
then if we introduce: 

D[p]=a\c-p), 
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we have 

_|C|2 1 

D[pm = = — (2.3) 

As a consequence, when a is large, a^(/9 — c) formany goes to A/? as for the non-local 
capillary term from [9], and the object of this article is to prove that the solutions of 
this system will go to the solutions of the local Korteweg model. Let us now define the 
interaction potential (/>q by: 

^(0 = ^^—^ (2.4) 
We have Jjj (j)a{x)dx = 1 and D[p\ = a'^{<j)a * P — p)- If we put </> = 01 then 

^(^) = TTiTTT' = ?(•/«)> and = a'^(l){a-). (2.5) 

In some cases we have explicit expressions for this inverse Fourier transform: for all x, 
4){x) = Ce~l^l when d = 1, 4>{x) = C'^-^ when d = 3 (we refer to [33]). In the other 
cases the expression of (f) involves Bessel functions. Let us begin by recalling that the 
fourier transform of a radial function is also radial, more precisely (see for example [33] 
page 213) there exists a constant Cd such that if /(x) = /o(|2;|) for all x G W^, then its 
Fourier transform satisfies for all ^ E M"^, f{0 = -^o(lCl) where for all p > 

^o(p) = / Jd^i{pr)h{r)r2dr, 

p2 ^ Jo 2 

where J^^ denotes the general Bessel function of real index v. This formulation is related 
to the Hankel transform, we refer to the appendix for more details and properties on 
Bessel functions. Coming back to our problem, we then obtain that for all x € W^, 

<Pi^) = r J^-i{r\x\)-^dr. (2.6) 

|x| 2~-^ Jo 2 i + r 

And thanks to the Hankel-Nicholson integrals (we refer for example to [27] page 330 or 
[36] page 434), under the following assumptions: 

3 

a>0, Re{z)>0, -l< Re{iy) < 2Re{fi) + 



we have the identity: 



where denotes the modified Bessel function of the second kind and index v (also called 
Hankel, Schlafti or Weber function). This allows us to finally write that for all x £ M*^ 
provided that d E {1, 2, 3, 4} (from the previous conditions with i/ = ^ — 1, /i = 0, z = 1, 

a = |x|), 

m = -^i^.„i(|x|). (2.7) 



X 2 
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Remark 3 Another way to understand the hmitation on the dimension consists in ob- 
serving in the integral (2.6), that if we roughly approximate the Bessel function by 
cos(r)r~^/^ at infinity, then the integrated function has the following asymptotic expan- 
sion at inhnity: cos(r)r'^/^~^/^ (\x\ = 1 for more simplicity). 

In fact (2.7) is also valid for dimensions d > 5 (Like the Fourier transform, the Hankel 
transform can be generalized for tempered distributions). Let us compute the Fourier 
transform: for all ^ € M"^ 

Kd_.{\x\) r°° d f r \ 

"""^^-d dx= / r2Kd_-^{r)( / e'^'^'^-^da; dr, 

and thanks to the radial symmetry: 

Performing a d-dimensional spherical change of variable we obtain that (with A = rl^l): 

-iAwei _ / / / ^-*Acos6»i • d-2 /i^ ^- c(-3 , 



-1 



JO JO 



e-'^'^-'^duj = / ... / e-'^'°'"'sm''-'eism^-'e2...smed^2d9i...ded^i 



Cd r e^*^™"^i sin'^"^ 
Jo 



Using the following integral representation of function ly (for Re{u) > —1/2 see the 
appendix for modified Bessel functions Ii, and Ki,) 



Iy{z) = ^ / e-^^^^^sin-^-^tdt. 

2^7r2r(z^ + i) Jo 



2 

Then thanks to the following identity (here a = and 6=1): 

zIy(az)Ky(hz)dz = -r (aIu+iiaz)Ku(bz) + bIi,(az)Ki,+i(bz)) , 

— 0^ 

we obtain that: 

. Kd_^i\x\) 2i-ir(|) 
^4 dx 



so that in (2.7), Cd = (22~^r(|))~^, Considering the asymptotics of function Kd_-^^, 4> is 

continuous on — {0}. Near 0, and for d > 3, we have 4'{x) ~ Cd\x\'^~'^ so that |x|(/>(x) 
is a function on W^. 

2.2 Reformulation of the system 

We are now able to write the system into a non-local form: 
{NSOPo) 

J dtPa + div {paUa) = 0, 

I dtipaUa) + div {paU U^) - AUa + V{P{pa)) = PaKa'^^{(f>a * Pa - Pa), 
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with 

(j)^ = a'^(t){a-) with (t){x) = — ^Kd^^{\x\). 
Remark 4 From the previous computations, we immediately get that 

Ca- Pa = (-A + ald^^Apa = (t>a * Pa - Pa 

that is Ca = (pa* Pa- This is why we cannot choose any initial data for the order parameter 
and take cq = (j)a * Po- 

As we consider initial data close to an equilibrium state (p, 0) we begin with the 
classical change of function p = p(l + q). For simplicity we take J) = 1. The previous 
system becomes (also denoted by (NSOPa))- 

(NSOPa) 

j dtqa + Ua-Vqa + (1 + qa)div Ua = 0, 

[ dtUa + Ua-VUa - AUa + P'{l).Vqa - Ka'^V {4>a * qa - qa) = K{qa)-'^qa " I{qa)AUa, 

where K and / are the real-valued functions defined on M given by: 



K{q)=(^P'{l)-^^^^^ and /(g) 



q+l 

The functional spaces we will really use are the following: 
Definition 2 Tiie space is the set of functions (q, u) in 

endowed with the norm \\{q,u)\\Eg = \\iq,u)\\E^{oo) where for all t we denote (recall that 
uq = min(//, v)) 

+ MMlib^+^' + ^h\\L\B'+'-^~^ + A\q.\\L\B'+^- (2-9) 

Remark 5 Due to obvious simplifications we slightly changed the notations for and 
-Bq^^'*: with the notations from [9] these spaces would have been respectively denoted 
byE' andB'j^''. 

l/a l/a 

We will now follow the tracks of [7] and [9] to prove the results. Classically in the study 
in critical spaces of compressible Navier-Stokes-type systems (see [12, 6, 18]), the proofs 
of theorems 2 and 3 (see [7] section 2) rely on key a priori estimates on the following 
advected linear system (a > is fixed and for more simplicity we write (g, u) instead of 
{qa,Ua)): 



(LOPa 



dtq + v.Vq + div u = F, 

dfU + v.Vu — Au + pVq — KcP'V{(j)a * q — q) = G. 
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With 

Au = ij,Au + (A + ^) Vdiv u. 

Although the potential function is different from the gaussian from [7], we can easily 
adapt the energy methods and results from this paper. Here we will directly focus on 
more refined estimates as in [9] and use them in the proof of the last theorem: we can 
prove that the estimates are similar up to slight changes in the constants: 

Theorem 4 Lei a > 0, -f + 1 < s < f + 1, / = [0, r[ or [0, +cx)[ and v G Li(/,^|^^)n 

d 

L^{I,B2i). Assume that {q,u) is a solution of System (LOPa) defined on I. There exists 
oq > 0, a constant C > depending on d, s such that if a > ao, for all t £ I (denoting 
u = fj, + 2X and uq = min(i/, /x) ), 

C _.C,,sc [\\\Vv{T)\\d +\HT)\\\)dT 

where 

1 4k; 
C 2 = Cmax(VP, — )max(-2, ( — )^), 

Cmsc = + l-^ + ^l +/^ + ^ _^ max(l, ■\). 

Remark 6 The coefficient Cmsc satisfies: 

_/i±^+max(l,i,) JfA + /i>0, 
li±^+max(l,J,) JfA + ^<0 

and when both viscosities are small, we simply have Cmsc < max(l, ^). 
2.3 Linear estimates 

As in [6] and [9] the first step to prove theorem 4 is to obtain estimates for the following 
linearized system: 



(OPa) 

With 



dtq + div u = F, 

dtu — Au + pVq — Ka^V{(j)oi * q — q) = G. 



Au = fiAu + (A + ^)Vdivu. 
In this article, we will use the following frequency-localized estimate: 
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Proposition 1 Let a > 0, s £ M, I = [0,T[ or [0, +oo[. Assume that {q,u) is a solution 
of System (Os) defined on I. There exists oq > 0, a constant C > depending on d, s, 
Co and Cq such that if ct > ao, for all t £ I (we recall that uq and C ^2 are defined in 

the previous theorem), and for all j € Z, 

\\Aju\\LocL2 + iyo2^^\Aju\\^L2 + (1 + i^2->) (^\\Ajq\\Loo^2 + 1^ mm{a'^ ,2^^)\\Ajq\\^L2^ 
< C^,^ ((1 + ^^2^')l|A,-go||L2 + \\A,uo\\l2 + (1 + u2^)\\AjF\\lil^ + \\AjG\\lil^) (2.11) 



(2.12) 



Remark 7 Let us precise that these linear estimates are valid for all dimension, the 
limitation d < 4 only appears in the advected case. 

2.3.1 Eigenvalues and eigenvectors 

In this article, since the methods are very close to [7] and [9] we will only point out 
what is different and refer to these articles for details. As in [12] or [6] we first introduce 
the Helmholtz decomposition of u. Defining the pseudo-differential operator A by A/ = 
^-i(|.|/(.)), weset: 

{V = A^-'^divti, 
w = A^-'^curln 

then u = — A~^Vf + A~^divt(; and the system turns into: 
' dtq + Av = F, 

(L'q,) < dtv - vAv - pAq + Ka'^A{(j)o, * q - q) = A"Miv G, 

_ dtw — fiAw = A~^curlG. 

The last equation is a decoupled heat equation, easily estimated in Besov spaces (see 
[2] chapter 2). Moreover, as the external forces appear through homogeneous pseudo- 
differential operators of degree zero, we can compute the estimates in the case F = G = 
and deduce the general case from the Duhamel formula. So we can focus on the first 
two lines and compute the eigenvalues and eigenvectors of the matrix associated to the 
Fourier transform of the system: 

The discriminant of the characteristic polynomial of A[^) is: 




A(e) = ierkier-4(p+« 



+ 1 

and thanks to the variations of function 



X , 2 A/ 2\ ^I^Op' 



fa'-x^vx — 4(p + K— ) = V X — A{p + Ka ) + 



^ + 1 ' 4 + 1 ' 
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we obtain the existence of a unique threshold Xq > such that 



<Oif <X„, 

> if > x„. 



We emphasize that this function has the same variations as in the case of [9]: when 

2 2 

3^ > 1) /a is an increasing function on M+, and when < 1, /„ is decreasing in 
[0,a^(^^^ — 1)] and then increasing. 

Proposition 2 Under the same assumptions, we have: 

4p 



a— >oo 



z/2 - 4k 



4k ' 



Proof: It is simpler than in [9] because here we have explicit expressions for the 
threshold: if we put A = u'^ — k — Ap/a^, then 



a 



2i/2 



A+\ 16p^ +A] . 



Next introducing the following function, we obtain the expressions of q and v exactly as 
in [9]: 

fa{x) 



9a{x) 



4 , X , 

z/^x ^ + 1 



(2.13) 



-For the low frequencies (A < 0), when |^| < yfx^, we have: 

m = ^ ((1 + ^y^- + (1 - ^Y'^-) m) - 

with: 



and 



4 |^|2 



+ 1 



(2.14) 



ii±iSiO). 



-For the high frequencies (A > 0), when |,^| > y'x^, we have: 

m = ^ ((1 - T^y)^*'" + (1 + T^y)^*'-) *^(^) + '^^TO^^'^o(^)' 
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with: 



m) = VaJm = - + ^^^) (2-15) 



and 

A± = -^(l±i?(0)- 
Remark 8 As in [9] it is crucial for the time integration to observe that 



2.3.2 Thresholds 



As in [9] we can find another threshold frequency Ua > Xa of size (in each case for 

2 

j^) that will enable us to push the parabolic regularization until frequencies of size a. 
In the present paper we will have explicit expressions. Let us first remark that rewriting 
function ga into the following form immediately implies that this is an increasing function 
from [0, oo[ to ] — oo, 1[: 



ga[x) = l 2~~T~T — \ — 2 withM= — 



(2.16) 



If /3 G [0, 1[ we easily compute that there is a unique positive solution of the equation 
Qaix) = f3 given by: 



1 



-A + 



16p 



SO that we immediately have the following asymptotics when a is large: 

— r-l|a2 ifM<-^, 

M(l-/3) J 1-/5' 



Xa,l3 ~ < 



P 



a 



iy\j l-l3 
ip M 



if M 



1 

1-/3' 



ifM>^. 



Remark 9 The previous proposition is obviously a particular case of this result. 
We are now able to define the second threshold Ua- 



-ifM<l, 2 
2 ~ v 

where M = — . 

l-^>lifM>l, ^'^ 
2M - 2 - ' 



(2.17) 



Using the previous result for /3 = ^ or 1 — according to the case for M, we obtain 
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Proposition 3 With the same notations we have that: 



If M > 1, Ua ~ a and for all a, 



<ya< 2a^ 



If M <l, Ua ~ (id- - l)a and for all a, 



M 2' 



2.3.3 Point wise estimates 



Now that we have defined the frequency thresholds Xa and Ua we have the following 
estimates. Up to the values of m and the second exponential from the density in the 
second case, they are the same as in [9] to where we refer for details or proofs: 

Proposition 4 Under the previous notations, there exists a constant C, such that for 
aU j G Z and all ^ G 2^C where C is the annulus € ]R'',co = | < |C| < Co = §}, we 
have the following estimates (we denote by fj = Ajf and we refer to the appendix for 
details on the Littlewood-Paley theory): 



• If 1^1 < ^: 

(1 + iy2^)mo\ < ce-^ ((1 + ^2^)1^(01 + (1 + jp)\^m) , 



[ m)\ < Ce-^ ((1 + u2^){l + VP)(1 + + 15^(01) ■ 

If^/x^ < 1^1 < y/y^: 

' (1 + vT)\qom < T^e-^(i-'") ((1 + ^2^)1^(61 + (1 + ^)l^(C)l) 



[\v,{0\ < T^e-^(i— ) (^.2^15^(01 + 1^^(61) 



wLere m = ^gaiVa) = ^ if M = <l, m = - ^ if M >l. 

(1 + i^TMm < C [e'"-^ + e-i:"^*) ((1 + z.2^-)|^(e)| + (1 + ^)l^^(C)l) , 

< C (e-^ + (1 - y^:(^)e-^(i-^/^^(^)) (i^2^|g^(OI + 1^^ 
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2.3.4 Time estimates 

As in [9], due to the choice cq = 3/4 and Cq = 8/3 (see the appendix), we can observe 
that there exist at most two indices = j^, — 1 or = such that E 2^[co, Cq] 
for j G ii^Ja}- 

We refer to [9] for the proof of the following proposition that implies Proposition 1. 

Proposition 5 Under the same assumptions as in Proposition 1, there exists a constant 
C such that for aU j eZ (denoting M = ^): 

• For aU j < j„, 



IvjWl^l^ + i^2^^\vj\\LiL^ + (1 + 1^2^) (Wqjh^L^ + '^2^'\\qj\\LlL^) < 



Cmax( — ,M^) (^(1 + i.2^)(l + ^)\\qojL2 + (1 + ^)\\vo,j\\L^ 



For aU j > j„, 

Cmax(l,M) (^(1 + u2^)\\qojL^ + (1 + -^)\\voJl^ 



2.4 Advected linear estimates 



(2.18) 



(2.19) 



The difficulties and methods exposed here are the same as in [9], so we will roughly 
explain them and focus on what is new. 

In order to prove Theorem 4, a natural idea is to use Proposition 1 and put the 
advection terms as external forces. Unfortunately, there are some obstacles: the main 
problem is that in v-Vq, the term T^^q can be estimated in -62^^ but not in i?| ^ because 
it is not enough regular in high frequencies. 

A direct use of the linear estimates will be useful only for the low frequencies [j < 0), 
and in the high frequency regime {j > 0), we will perform a Langrangian change of 
variable (as in [22], [23], [15], [6], [9]) in order to get rid of v • Vq. We then aim to 
use on the new system our linear estimates but we have to be careful with the external 
force terms introduced by the change of variable. Most of the work in [9] was to provide 
estimates on the commutator of the non-local operator from the capillarity term and the 
Lagrangian change of variable. 

For all j € Z and t £ I, we introduce: 

Uj{t) = \\Aju\\L^L^+uo2^^\\Aju\\^L2 + {^+u2^) (^\\Ajq\\L^L^ + umin{^,2'^^)\\Ajq\\LiL^^ 

(2.20) 
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and 



Uit) 



2,1 



LfB. 



2,1 



(2.21) 



2.4.1 Low frequencies 

For the low frequencies, we obtain (see [9] section 3.1 for details) that there exists a 
nonnegative summable sequence whose sum is 1, denoted by (cj(t))jgz such that for all 
J < 0, -ftr > (to be chosen later), and if a > A'^i/log2 (we refer to [9] section 3.1 
for this, and to the appendix for A'^i which is a constant related to cq and Cq in the 
Littlewood-Paley decomposition, such that if |j — Z| > A'^i then Aj o = 0, in our choice 
of Co and Cq, A^i = 1). 



C/,(t)<C 



p.- 



?7,(0) + (1 + z.2^-)||A,-F||^,^. + ||A,G||^ii. 



'^O + l + l^||(7||^a + l,s-l -\- V ||(/||^s + 2,s 



2 Jo 



Cj{T) 



■max(l,^) + — )||7;(T)f , + C\\v{t 



dT 




Remark 10 As pointed out in [9] (remark 29) using the linear estimates for the low 
frequency case, allows us to get rid of another difEculty introduced by the change of 
variables: in low frequencies some of the additional external force terms have too much 
regularity to be absorbed by the left-hand side, and not enough regularity to be controlled 
with a view to apply the Gronwall lemma. The only way to control them would be to 
use interpolation arguments that would introduce linear time dependant coefficients, and 
prevent us to get global in time results. 



2.4.2 Lagrangian change of coordinates 

As explained, in order to get rid of the advection terms involved in system {LOPq) the 
first step is to consider the following localized equations (as usual we set fj = Ajf...): 

dtQj + Sj-iv.Vqj + divuj = fj, 

dtUj + Sj-iv.Vuj — Auj + pVqj — Ka'^V{(l)a * qj — qj) = gj, 

where the external force terms are defined by: 

fj = Fj + (^Sj^iv.Vqj — Aj{v.'Vq)^ and gj = Gj + (^Sj^iv.Vuj — Aj{v.'Vu)^ . 

Both of these terms can be estimated thanks to the following commutator estimate from 
[15] (we refer to lemma B.l from appendix B): 
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Lemma 1 ([15]) There exists a sequence {cj)j^z £ /^(Z) such that ||c||;i(2) = 1 and a 
constant C = C{d, a) such that for all j E Z, 

\\Sj-iv.Vhj - Aj{v.Vh)\\L2 < Ccj2-^''\\Vv\\ . d \\h\\^. 



In order to perform the change of variable we define ijjj^t as the flow associated to Sj-iv: 

(dttpj,t{x) = Sj-iv{t, ^2 23) 

we can also write: 

ft 







Thanks to propositions 8 and 9 from the appendix (we refer to [15] or [6]), there exists 
a constant C such that: 



US' o V'i.t ^ c'^^II^IIlp for all function g in L^, 
m%-Id\\L^ < e^^-l, 

ll-D'^V'tji- < C2('=-i)j(e^^ - 1) for A; > 2, 



(2.24) 



where 



V{t) =^ [' \\Vv{T)\\L^dT. (2.25) 







As explained in [9], considering some of the additionnal external force terms at the 
point '4}~l{x) instead of x will help getting uniform estimates with respect to a, so the 
jacobian determinant of the change of variable will play a crucial role in our computations 
(contrary to the case of lemma 2.6 from [2] where it produces a term that we are not 
able to sum here). 



where Xj{T,t,x)) denotes the two parameter flow associated to Sj-iv (we refer to (4.88) 
in the appendix). 

Let us now perform the lagrangian change of variable, for a function h, we define 
h = ho ipj^t = h(t,ipj^t)- Then we have dtqj{t,x) = {dtqj + Sj-iv.'Vqj){t,^j^t{x)), which 
provides the following system: 



dtQj + diviij = fj + Rj, 

dtUj — Auj + pVqj — Ka''V{(f>a * Qj — Qj) = dj + Rj + Rj + kRj, 
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where the remainder terms Rg, Rg and R^ are exactly the same as in [6] and [9] (with 
the same convention: if / : M"^ — t- is a differentiable function then Df denotes the 
Jacobian matrix of /, and V/ is the transposed matrix of Df.): 

R]{t,x) := Tr{Vuj{t,x) ■ - Vi^Tli^P,,tix)))) , 

R]{t,x) := Vqj{t,x) ■ {Vi^jjii^jAx)) " ^d) 

and R^j := fj.R'^ + (A + i^l)R^j with 

Rf\t,x) := Tr((V^7,^(^j-t(x)) - h) ■ VDu%x) ■ Di^-jiiPj^x)) 

+VDu]{t,x) ■ {DiprJiiPjAx)) - Id)) + Vu%x) ■ A^-/(^,- t(x)) 

R]'\t,x) := Tv{Duj{t,x) ■ duD^^^i^j^x))^ 

+ dl^u]{t, x) ■ dkTp~l'''{i^j,t{x)) ■ daip'j'^iipj^x)) 

a,b,c,b^a,c^k 

d 

+ 'Ydliu){t,x) ■ {dkip~}'^{t,^IJj,t{x)) - Id) ■ diipjl'\t,'ipjAx)) + {diipjl^'i^j^x)) - la) 

i=l 

As in [9], the only difference with [6] is the following additionnal remainder term: 

Rj = a^{(Pa * Vqj - Vqj) o V'j-j - a^((^„ * Vqj - Vqj). (2.28) 
Thanks to the definition of (pa and (j) we can also write that: 

LM) '= o^H4>a *f-f) = o? f m (fix --)- fix)) dz. (2.29) 

As in [9] most of the work consists in obtaining bounds in L^L^ that are uniform with 
respect to a, and go to zero when t is small and dealing with Rj is the object of the rest 
of this section. 

2.4.3 Precisions on the capillary term 

Let us first go back to the convolution term written in (2.29): for a function /, 

a'(<^.*/-/)(0 = -^^/(e). 

Similarly to [9] we obtain the following equivalence, giving a smooth interpretation of 
the hybrid norm. Here again, instead of a fixed frequency threshold there is a continu- 
ous transition from the parabolically regularized low frequencies and the damped high 
frequencies: 

Proposition 6 For any suitable function f and any s gM, we have: 
.+2,. = J]min(a2,22^-)2^-^||A,/||i2 
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Proof: Thanks to the monotonicity of function x i— >■ x/(l + x), we easily obtain that for 
all j G Z: 

1 2^"' 
We refer to [9] for details. 

Remark 11 In the U-setting, we can prove that for all j G Z: 

\\a\<Pe*^,f-^jf)\\Lr < Cmin(a2,22^')||A,/||i.. 



2.4.4 Estimates on the capillary remainder Rj 

This section is devoted to giving estimates on the capillary term introduced in (2.28). As 

Vqj = V{qj o ipj^t) = Vqj o iPj -t x Dijjj^t = Vg^ o * x {Dipj^t - Id) + Vg^ o i/jj^t, 
we obtain the following decomposition: Rj = Ij + Ilj with 

Ij = a^icpe * gj - gj) where gj = Vqj o ipj^t x {Id - Dipj^t), 
Ilj = a\(l)e * Vijj - Vqj) o ijjj^^ _ f * {Vqj o - Vqj o Tpj^t j . 



So that as in [9] , we need to estimate (locally in frequency) the commutator between the 
Lagrangian change of variable and the non-local operator Lq defined in (2.29). For a 
function /, and all j G Z we set fj = Ajf and: 

Il'j = Il'j (/) = (0, * -fj)o i^j^t - * (/, o ijj^,) - fj o Vj,t) . (2.32) 

Theorem 5 Let o" G M. There exists a constant C = C(j4 such that for all f G Ba^^''^ , 
there exists a summable positive sequence {cj{f))j^z whose sum is 1 such that for all t 
so small that 

e^^^ -l<\- (2.33) 

and for all j G Z, 

||/7;.(/)||^. < Ce^^(V + e^^^ - l)c,(/)2-^ni«'(</'a * / " 
where V{t) = \\Vv{T)\\Lo^dT. 

Remark 12 As a by-product we obtain that under the previous assumptions, if t is 
small enough, 

a^{<|>,*Ajf)o^p^,,-cP,*{Ajfoi;j^,)\\^,<Ce^y{V+e^<'''-l)cj{f)2-^^ 
Remark that neither of the left-hand side terms are spectrally localized. 
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Proof: we refer to [9] for details and here we will only focus on what changes. The first 
step is to obtain pointwise estimates and then LP' estimates: as in the works of T. Hmidi, 
S. Keraani, H. Abidi and M. Zerguine ([22], [23], [24] and [25]) we wish to retrieve the 
desired Besov norm thanks to an equivalent expression of 11'^ as an integral formulation 
involving finite differences of / of order 1 (that is expressions of the type T-yf — f where 
T-yf{x) = f{x + y)) or order 2. Like in [9] we need to directly consider I I'-lil^jl {x)) 
instead of simply IIj{x) and we obtain: 

Il'.i^Pjjix)) = a' [ <P{z) (fj{x - ^) - fj{x)) X 

0(a(V,7l(x)-V^7l(x-f))) - f idivS,^iv){T,XjiT,t,X--))dT 



dz. (2.34) 



Remark 13 We emphasize here that the previous quotient is weU defined near zero and 
we refer to section (2.1). 

We have to face the same problem as in [9]: as we want to estimate the norm of this 
quantity, that is a Besov norm with integer regularity index s = 0, the finite difference of 
order 1 will not be sufficient for our need, and we will have to introduce finite differences 
of order 2 (this is a classical problem for integer indices). Indeed, using the present 
quantity would only involve a term in \y\/a and when estimating in low frequencies, 
there would be either an extra multiplicative coefficient a or an extra derivative term 
(that would prevent any convergence when —j is large). To be able to do a correct 
estimate we need at least |?/p/a^. 

For this, we simply write II j = \{II'j + ^I'j)^ ^-^^d perform the change of variable 
z = —y in the second integral. If we set: 



B. 



0(a(^7i(x)-^-i(x-f))) 



0(a(V^-i(x)-V'-i^(x + f))) 



(2.35) 



z , 



(div5j_iv)(T, Xj(r, t, X ))dT, 

" 



{div Sj-iv){T,Xj{T,t,x H — ))dT, 
" 



then 



Il',{^-}{x)) = ^ / <t>{z) (Mx --)- fj{x)) [1 - R^e-^-]dz 



+ 



a 



a 



fj{x) ) [1 - R+e''^+]dz = IIIj{x) + IVj{x). (2.36) 
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where 

„2 



III,{x) = ^ / 0(z) (fj{x --) + f,{x + -)- 2/,(x)) [1 - i?„e-^-]d^, 

2 (2-37) 

I 2 Jjjd V a / 

Taking L? norms we have: 

II//;. o < \\iiij\\L2 + ii/vs-iil2, (2.38) 

and thanks to estimates on the Jacobian determinant of the flow (see 2.24), theorem 5 
is immediately implied by the following proposition: 

Proposition 7 Under the previous assumptions, there exist a positive constant C = C„^d 
and a nonnegative sequence {cj = Cj{f))j^z whose sum is 1, such that if t is so smaU 
that e^^^W - 1< i we have: 



\IIIj\\L2 + \\IV,\\l2 < C{V + e^^^ - l)e^^2-^''c,\\a\cl)a * f - f)\\ 



^2,1 



To prove this result we will successively prove the following lemmas: 

Lemma 2 There exists a constant C only depending on the dimension d such that for 
aU j G Z, /, and all t is so small that e^'^^^*) - 1 < ^, 



|///,||i2 <C7e^^(y + e^^^ -1) 

xa^/ e^^ ^ ^,,3 ' m- --) + M- + -) - 2/i(.)llL^^^, (2.39) 



and 

||/y,||i2 <Ce^^(y + e2^^-l) 

1 • 

xa^/ e-^ ^ + ' min(l, ^)||/,-(. + -) - f.^L^dz, (2.40) 
where V is defined in (2.25). 

Lemma 3 For all a & M there exists a constant Ca,d such that for any f G B21, there 
exists a nonnegative summable sequence (cj(/))jgz with ||cj(/)||;i(2) = 1, such that 

I ^~'^^^Tr^\\M---) + f^(- + -)-^fM\L2dz 

jRd \z\'^-2 ' a ' a 

I I — i '11 

+ " / ^ ' , ,^3 mm(l,^)||/j(. + -)-/j(.)||^2rfz 

< a,rf2-^%(/)||a2(<^, * / - /)||^^.^. (2.41) 
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Proof of lemma 2: A direct estimate gives: 

2 f 

< ^ / 4>{z)\\f,{- --) + fji- + -)- 2/,(-)||l2 111 - R.e-''-\\L^dz, 



and 

,2 



\\IV,\\l^<'^ / 4>{z)\m- + -)-fmLl\\R-e~''--R+e-^nL^dz, 
So that we first have to focus on the L°° norms: 

Lemma 4 Under the same assumptions, for all rj €]0, 1[, there exists three constants C, 
C(i and Cfi^n such that for all j G Z, all f, and all t > so small that e^^'^i^) _ i < 1^ 
have: 

111 - iJ-e-^iUgo < C,,,e^'^^W(e2^^ - 1 + V)e\'\\^+^^-^t'"'''''''-%l + |z|-|). 



a 



< C,,,e^^^W(e2^^ -1 + V) min(l, ^)el^lK(i-^)(^^^"^^^-i)] (1 + |^|-|). (2.42) 



Proof of Lemma 4: as in [9] every term is close to 1 when t is small and we obtain the 
result by carefully estimating the differences between these terms. For this we write: 

\Q2 = R-e^^- -R+e~^+ = {R.-R+)e"^' +R+{e~^~ -e^^+). ^ ' 

Thanks to the following elementary consequence of the mean-value theorem: 
Lemma 5 For any G M, |e^ — e^| <\x — y\e^^^'^^'y\ 

we obtain (we refer to [9] for details) that there exists a constant C > such that (we 
refer to theorem 5 for the definition of V): 



\l-e-^-\<CV{t)e''''^'\ 

\z\1K 



le^^-Ke^^^*)^ (2.44) 



-i?- _e-^+| <2Cy(t)e2^^W min(l, 



a 



We now turn to the estimates involving R± . As obtained in (2.7) , the interaction potential 

can be expressed thanks to the modified Bessel function of second kind Kd_^ : 

2 ^ 

'^(^) = --§-r-^i^-i(kl)' 

X 2 2 



SO that if we denote by Y± the variation ratio of the flow: 
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we have the expression: 



2 2 

We bound the first factor thanks to the following estimates on the variation ratio: 



(2.46) 



Lemma 6 ([9] section 3.4) Under the general previous assumptions, for all x,z £ 
with z ^ 0, we have: 

^-cv{t) ^Y^< e^^W 



and 



_ 1| < e^CVW - 1, 



Y± 



< e^^^W - 1. 



Using the upper and lower bounds for cj) given by proposition 11, we can write that for a 
fixed r] g]0, 1[ (to be precised later) there exists a constant Cd,r] such that: 



C(l-l)V{t)_ 



{l-rJ)a\^(,-lix)-^(,-lix±^) 



R± < Cd,r^e 

(a|V,-/(x)-V-i'(x±f) 
that is, with the previous notations: 



d-l 
2 



|2;|2-^ ifd>2, 
{l + \z\^) ifd = 2. 



R± < Cd,r,e 



C{l~l)V{t)lz\il-{l-v)Y±) 



d-l 

1 \~ 1 



Y± 



\Z\ 2 



if d > 2, 
if d = 2. 



and then 



R± < C,,,e^('^-i)^We^l^le(^-^)l^l(^-^±)(l + |z|-|). 
Thanks again to lemma 6, we finally obtain: 

R± < Q,,e^''^Wel^lK(i-'')(^'^"*'^-i)](l + |z|-i). 



(2.47) 



There two terms left to estimate: 1 — R- and i?„ — Again, we simply write that: 



1-R± = P1+P2 



1 \^ 



+ 



1 \^ 



Ka A\z\ 



An elementary use of the mean- value theorem to the function /i(y) = 2/2 ^ gives that 
(we recall that in lemma 6, we proved e"*-"^^*-* < ^ < 6*^^^*)): 



IMl) - Mt^)I < |1 - T^l(| - 1) sup <Cd\l-^ 



\^c\\-2\vit) ^ (2.48) 



that is: 



lAl <Q(e2^^W-l)e^''^W 



(2.49) 
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Moreover, the second term is bounded by: 



IP2I < 



oCdV{t) 



2 ^ 

Similarly, we can write that: 

Ki^,{\z\) - K,_,ia\^-l{x) - ip-J{x ± -)|) 



K.,{\z\) - Ka,{a\^-j{x) - ^p-lix ± ^)|) 



<|z||y±-l| / K',^^{\z\{l + u{\Y^-l\)) 



du 



y<^ [l^l (i_(e2Cy{t)_i)) (i+(e2CV(t)_i)) j 



. (2.50) 



If t is so small that e — 1 < 2' then thanks to the bound for K^^ from proposition 

2 ~ 

11 we obtain that: 



K,_^{\z\)-K,_^{a\^jl{x)-i;;l{x±-)\) < (e^^^W - l)|z| • sup Ca, 



?/6[i|z|,f|z| 



y 2 



(l_^)|^|(l_(e2CV(t)_i)) 

<C7,,,(e2^^«-l)^ . (2.51) 



\Z\ 2 



Gathering this with the lower bound for Kd_, from proposition 11, we finally bound P2'- 



IP2I < e^''^WQ.(e2^^W - 1)- 



-(l-„)|.|(l-(e2CV(t)_i)) 



\Z\ 2 



(1 + kl 



if d > 2, 
if d = 2 



< Q,,e^'^^W(e2^^W - l)e'?l-le(i-^)l^l(^''^'''"-i)(l + |z|-|). (2.52) 
Finally together with (2.49), we finally obtain the estimate: 

|1 - R±\ < C,,,e^''^W(e2^^W - l)el-lh+(i"'?)('^^^^<^^-i)](l + |z|-|). (2.53) 
Let us finally turn to the last term: i?_ — R^. Using the notations, we can write: 
1 X f-i Kd,{\z\Y^) / 1 \ f-i Kd,{\z\Y+) 



R.-R4-=\ — 



1 



Y: 



Ki_,{\z\) 



1 \ 2 1 
1^ 



i^l-i(l^l) 

= P{+P^. (2.54) 



As we did in (2.48), if k{y) = y we can write: 

|M^)-fc(^)|<Q|y_-y+|e^''^w, 
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The other term in P[ has aheady been estimated when estimating R±: 

Kd_^[\z\) 

2 

and we get that: 

\Pi\ < Q,,e^^^W|y_ - y+|el^lMi-'')(^^^"'*'-i)](l + Izj-I) 
For P2, as in (2.50), we can write: 

K,_,{\z\Y.) - K._^{\z\Y+)\ < \z\\Y^ y+l K',_^ (^|z|((l - u)Y+ + uY^)^ 
< \z\\Y. - y+l 1^ K',_^ (^\z\ + \z\{{l - u){Y+ - 1) + u{Y^ - 1))) 



(2.55) 



du 



du 



<|z||y_-y+|- sup 

j/6 [l^l (l-(e2CVW_i)) ,1^1 (i+(e2CV(t)_i)) j 
<Cd,r,\Y.-Y+\ 



d-l 
\Z\ 2 



(2.56) 



Using once again the lower bound for Kd_, imphes that: 

2 ^ 

l^'^l < C,,,e^^^W|y_ - Y+\e\^\V>+^^-^^^-'^"'''-^)Xl + \z\~h 
and together with (2.55), this implies that: 

- R+\ < Q,,e^'^^W|y_ - y+|el-lh+(i-'?)(<=^^"'*^-i)](l + Iz^l). 
Thanks to lemma 6, we get: 

|y_-y+|<e^^W(e2^^W-l). 



(2.57) 



(2.58) 



Unfortunately, as explained, this will not be sufficient: this estimate is useful for high 
frequencies j, but after integration in z, the result is not summable when j goes to —00. 
This is why, as in [9], we need a much more precise estimate on |y_ — y+|: 



|y_ - y+ 



a k,-/(^) - - f )!' - \r,l{x) - jjjjjx + f )P _ g Nj 

\z\ |^7i(^) - ^Tl(x - ^)| + |V;-i(x) - V-/(x + f )| \z\ D, 



(2.59) 



As in the proof of lemma 8 from [9], we simply use the identity |ap — |6p = {a + b\a — b) 
and write: 

iV, = (2^7i(.) - ^-i(. -1)- ^-i(. + -1)- ^-i(. + ^)) 

(-) - ^Jii- -"-)- ^jK- + ^)I • - ^) - + (2-60) 



< 
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Thanks to the mean value theorem (used twice for the first factor and once for the 
second), we can write: 

iVj < {—?\\D%1\\l^ X {\^)\\Di;-l\\L^, 
a ■'' a ■'' 

and using the estimates for the flow (see (2.24)), we obtain that 

iV, <e^^(e^^-l)^. 

On the other hand, 

n. = M(y^ + y_)>2e-^^wM. 
a a 

From these estimates, we easily conclude that: 

2^z\^ Inva.ocvuCV -2^l2l 



|y_ _ y^l < ^ X e^^ (e^^ - 1)^^ x -e"^^ f- < e'^^ (e^^ - 1) 

\z\ a-^ 2 \z\ a 

Then, combined with (2.58), we obtain: 

|y_ - y+l < e^^^ie^^^ - 1) min(l, ^). (2.61) 

a 

and then, thanks to (2.57), we obtain that: 

\R. - R+\ < C7,„e^''^W(e2^^ - 1) min(l, ffl)el-lh+(i-'')(-'^^'''-i)] (1 + l^^i). (2.62) 

a 

Gathering (2.44), (2.47), (2.53) and (2.62), we obtain the result from lemma 4. ■ 



End of the proof of lemma 2. Thanks to lemma 4, and proposition 11 (with the 
same rj), all that remains is to estimate the following function: 

cp{z)e\^\[nHi-v){e^^^'^^^-i)] < c^_^^:!^,N[.+{i-.)(e^^-(^)-i)] 

\z\'^ 2 

< ^^^^^_g|.|[-l+2,+ (l-,)(e^C.V(*)_i)] ^2.63) 
\z\'^ 2 

Using that e^'^^^*) - 1 < ^, we have 

-l + 2r? + (l-ry) {e'''''^'^-^)<-l + lv<-l 

if we choose rj = j [in fact we need < i). ■ 

Proof of lemma 3: as fj = Ajf, we can write that for all x and z, 

fj{x - ^) + fj{x + 8^) - 2fj{x) = (t_^A,/ + T^Ajf - 2A,/) (x) 

and we refer to [9] (proof of lemma 5) for the following result which is adapted from [2] 
(theorems 2.36 and 2.37): 

25 



Lemma 7 Under the same assumtptions, there exists a nonnegative summable sequence 
of summation 1, that we will also denote by {cj)j^z, such that: 



\T_^Ajf + T^Ajf-2A,f\\L2 

22j| -,12 

<a2^^'^Cjmax(Q-2,2-2j)||a2(<^„*/-/)||^. niin(l, (2.64) 

2,1 Q/'' 



< a2~^'^c, max(a-2,2"2j)||a2(<^^ * / _ f)\\^^^^ min(l, ^) (2.65) 
It appears then clearly that both integrals are bounded by: 

2 

a2-^%||a2(</.,*/-/)||^.^ max(a-^2-2^■)a2xI = C,2-^%||a2(<^^*/_/)|^ 
where we define: 

^def f _Ml + |z|-| . ,^ 22J>|2 p _.l + r-^ . ,^ 22ir2 

1= e 8 1— Lg — min(l, -^)dz = Cd / e s =— min(l, — 5— )?" ur 



Cd / e—s{l + r2)min{l,^^)dr = Cd(Ii+h), (2.66) 



where we have splitted the integral into: 

2^^' fij _r_ n, 1, f°° _r_, 1, 

Ii = — ^ / e sr (l + r2)(ir and I2 = / e 8(l+r2)(ir. 

As in [9], the second integral is easily bounded: as there exists a constant C > such 
that for all x > 0, (1 + r2 )e~8 < Ce~i6, we have 

e"T6dr < 16Ce"^^. 

2J 

For the second integral, if we use the same argument, we get the estimate: 

Ii < C7— / e-T?,dr < 16C^min(l, -r) (2.67) 
Jo 2-? 

which, multiplied by max(l, ^) gives a resulting term in max(l, ^) that is not summable 
for high frequencies. To overcome this difficulty we simply write that there is a constant 
C" > such that for all r e [0, §], 

o IrOi Ir (y. r 

r (1 + r2)e~8 < — r(l + r2)e~8 < C'—e~^. 

and then we obtain: 

2J Jo 
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which is obviously interesting only when ^ < 1 that is for high frequencies: in this case 
it cancels the diverging max(l, ^). But in low frequencies this new term is much bigger 
that the one from (2.67), so that we have to combine both estimates and finally get that: 



^22i ^ a, /"t _^ , ^22J ^ a2 

< (7^5- minfl, — ) / e i6 dr < 16(7^5- min(l, 
2-? Jo 



^1 



We are now able to estimate I: 



v2 



/2'^^ a" _J_JL 

I<C,(^^min(l,^) + e i6 

and then 

2 

C<,2"^-Cj||a2(</,^ * / _ max(l, ^) x I < 



22i 

_ • 9 ( o? 2^^ o? _ j_ ^ a" _ j_ ^ 

^2 J'^Cj||a'(</.„*/-/)||^. ( max(l,^)-^min(l,^) + max(e ^«2.,^e 2. ; 

< a2-J-Cj||a2(0„ * / - /)||^^.^, (2.68) 

which concludes the proof of lemma 3. ■ 
2.4.5 End of the proof of theorem 4 

As this part is strictly the same as in [9] (section 3.5) we will not give details (in particular 
we refer to [6] or [9] for estimates on ' and to [35]). Going back to system (2.27), we 
use the linear estimates from proposition 1: for all I G Z, as qj{0) = qj{0,'4>jfi{.)) = qj{0), 

W^lUjh^L^ + uo2'^^\\AiUj\\LiL2 + (1 + u2^) (^\\Aiqj\\L^L^ + z^min(^,22')||Ajgj||iii2 

(1 + u2')\\Aiqj{0)\\L2 + \\AiUj{0)\\L2 + (1 + i^2')\\Aj, + A/^H .1^2 



P' 4k 



+ \\Aigj + AiR] + AiR^j + KAiRjU^iL^ 



. (2.69) 



Thanks to (2.22), all we need is to estimate the high frequencies, that is {qj,Uj) for j > 0. 
For this, we define some A'^o G Z (that will be fixed later), and write: 

\ l>j-No 
We refer to [2] lemma 2.6, or [15], lemma ^.1) for the following classical estimates: 
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so that going back to Uj (we refer to (2.20) for the definition), we can write that for ah 

i>o, 



Uj{t) < Ce 



cv 



[e^v _ 1 ^ 2-^°e^^) Uj{t) + max(l, 22^«)C ^ : 



((1 + u2')\\Aiq,mL^ + W^i^mh^ + (1 + z^2')(||A;/,||^j^. + \\AiR]hiL, 

l>j~No 



(2.70) 



We refer to [9] for the estimates on the remainder terms R^^'^'^\ Recall that the previous 
section is needed to estimate Rj. If we fix A'^o > large enough and t so that: 



4 - 8 

^C(e^^-l) fl + 25^«C ^( 



(2.71) 



we obtain that for all j > 0, 



U,{t) < 3C25^°C ^ (1 + u2^)\\q,m\L^ + II%(0)||l2 



+ (l + z.2^-)||F,-||ij^. + ||G,-||iii.+ / 2-^^^~'^c,{t)\\Vv{t)\\ , U{T)dT 



Now, if t is so small that: 



3C25^«C ^(y + e2^^-l)<^-, 



(2.73) 



and if we take K = (2C ^) -"^ in (2.22), then sum over j £ Z, we end up with 



U{t) < ^ + C ( f/(0) + WFh.j^s-. + ^^Ili^lLiBi, + \\G\\ 



+ ( 



1 + |A + /i| + ^ + z^ 



1^0 



+ max(l,-^))^ Vt^'(r)[/(r)j (2.74) 



where 



F(t) f \\Vv{T)U^dT < W{t) [\\\Vv{t)\\ , + \\v{T)f , )dr. (2.75) 



^2,1 
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and thanks to the Gronwall lemma, we obtain that for t small enough (satisfying condi- 
tions (2.33), (2.71) and (2.73)), 



A + \X + n + v . 1 

Then we globalize the result as in [9] 



2C ,2{—^ — ^ +max(l,^))VF(t) 
xe ^0 y""" (2.76) 



2.5 Proof of theorem 2 

2.6 Existence and uniqueness 

As explained in the introduction, once we have defined the interaction potential (/>q , we 
can follow the very same methods as in [7]: using energy methods gives apriori estimates 
on the advected linear system, then the proof for existence and uniqueness is classical 
and follows the lines of [2] (section 10.2.3) for the compressible Navier-Stokes system 
(see also [14], [12]). In order to use the classical Friedrichs approximation, we define the 
frequency truncation operator J„ by: for all n G N and for all g G L^(M'^), 

JnQ = (l2-"<|C|<Co2"(05(6) > 

and the following approximated system (we omit the dependency in a for more simplic- 
ity): 

J dtQn + Jn {JnUn-'^ Jnqn) + Jndivu„ = 

I dtUn + Jn (JnUn-^JnUn) - AJnUn + P'{p).VJnqn " Ka^V(0a * J„g„ - J„g„) = G„, 

where 

\Fn = -Jn{Jnqn-div JnUn) 

\Gn = Jn{K{Jnqn)-^Jnqn " I {Jnqn)AJnUn) 

It is easy to check that it is an ordinary differential equation in x {L'^)'^, where 
= {n G L^(]R'^), JnU = u}. Getting uniform estimates implies global lifespan if the 
initial data is small enough. Then classical compactness arguments give existence of a 
global solution. 

Remark 14 Let us emphasize that the precise estimates proven in the present paper 
cannot be used in the proof of the existence: indeed the term Jn (wn-Vg„) is an obtacle for 
the lagrangian method to give bounded constants (with respect ton). So for the existence 
we simply use the rough apriori estimates given by the classical energy methods (as they 
rely on inner products in L^, here J„ has no effect). 

To obtain uniqueness for a fixed a, using (2.10) the computations are close to those for 
the compressible Navier-Stokes system. As in [2] we need to separate the cases d > 3 and 
d = 2 (the case d = 2 is more difficult because of endpoints for the remainder estimates 
in the Littlewood-Paley paradecomposition) . 
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2.6.1 Uniform estimates 

Using estimate (2.10) with s = ^, we obtain that for all t G M+ (we refer to [7] for the 
expressions of the external force terms F and G): 



g{t) = \\{qa,U^)\\ d <C 26 
E^{t) P'lT^ 



t 

2 



C ,,2C,isc / (||Vn„(r)|| , +\\u^{T)f a )dT 



Hi Hi 

X (/i(0) + , , + v\\F\\^ d +I|G|U ^-i) (2.77) 

^ L]Bl, \lBl, \}Bl, V 

Let 77 > be small (it will be fixed later) and assume that: 

5(0) =^ ||n(0)|| . + Um . + ^||g(0)|| . . < r?. 

-^2,1 ^2,1 -^2,1 

Let us now define 

r = sup{tGM+, <7(i) < 2C ^5(0)}. 

As g(0) < 77, we have T > (C > 1) and the aim is to prove by contradiction that 
T = 00. Assume that T < 00, then we obtain that for all t < T, denoting by C = C" 2 

and C^,jsc = Cvisc/ 

g{t) < Ce2C^^^»..^(i+2C.+e-^-i.c.(i+-.))^(0)_ 

2 

So that if = Ty(p, |-, Cmsc/^o) > is small enough then for all t < T, g{t) < 2C ^g{0) 

P'4i^ 

which contradicts the fact that T is maximal. It then implies that T = 00. 

Remark 15 When for example A = we have fi = u = uq and for a snaaU v, C'^j_g^ ~ 

so that the previous condition simply implies that r = rj/v^ < 1 must be small enough 

so that e2'C^^^(i+2C+e^'^ ) < 2. A sufRcient condition is that r satisfies e^^'^^^^'^ < 2. 

2.7 Order parameter estimates 

As already explained (see remark 4), we liave — (pet 

* Pa, SO that Ca — 1 = 4>a * Qa and 

for all s G M, IIcq, — IIIrs < IIooUrs . Moreover: 

-°2,1 -°2,1 

]^ 2^' 

\\Ca - PalU^ = * Qa - ^allns = ^||ga|Us+2,. = V'2'''min(l, ^) || Aig^ ||^2 . 

2.1 2.1 Q/-^ ' ^ (Y^ 

On one hand, thanks to the Lebesgue theorem (for series), for sG{^ — 1,^} this involves 
that: 

||Ca-/3a|U .4-1 +J^\\Ca-Pa\\-. d >0, 

and on the other hand, thanks again to the energy estimates, we end up with: 

2„ „ C° 

l^\\Ca-Pa\\ +l^\\Ca-pa\\ d 

Li{R+,B^-^ ) L1(R+,B|J Ct^ 

This concludes the proof of theorem 2. ■ 
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3 Rate of convergence (Theorem 3) 

In this section we prove that the solution of [NSOPa] goes to the solution of (K), and 
we give estimates of the rate of convergence as a goes to infinity. Here we follow what 
we did in [7]: everything relies on Theorem 4. As already explained, if the initial data 
satisfy 



|go|| .4-1 +i^lko|| d + \\uo\\ d_^<r] < mm{r]K,'nop), 

^2,1 -°2,1 



^2,1 



then systems (K) and (NSOPa) both have global solutions {q,u) and {qa,Ua), and with 
the same notations as before, denoting C = C ^2 then for all t G M we have, 



9a{t) =^ \\iqa,Ua)\\ d =^ \\u\ 
E2{t) 



<i_, + g L d_i+i^ k L d 

L^BL L^Bl, L^Bl^ 



+ . 4+1 + T^h\\~, . 4+1,4-1 + ^'^IklU, . 4+2,4 < Cr/, (3.78) 



^2,1 



and 



ffHi) = ||n| 



d-i + \\q\\~ . 4-1 + 1^ 

r 00 D 2 



-'^t -°2,1 ^t-°2,l 

+ . ''+1 + ^^Ikll .d+2^C!r7. (3.79) 



As in [7], up to an additional forcing term, let us rewrite system [K) with a capillary 
term as in system (NSOa)'- 



dtq + u.Vq + (1 + g)div n = 0, 

dtu + li.Vn -Au + P\l).Vq - Ka'^V{(j)a * q - q) = K{q).Vq - I{q)Au + Ra 



where the remainder Ra kV (Ag — a^(</>Q * q — 5)) and K and / are defined in the 
introduction. Let us now write the system satisfied by the difference {6q, 6u) = {qa — 
q,Ua - u): 



dtSq + Ua-^Sq + div 5u = 5F, 

dt6u + Ua.V5u - A5u + P'{l).S/5q - Ka^V{ 



6q) = 6G-Ra, 



(3.80) 



where 



5G 



dg 

def 



Ell SG, 



with < 



6F3 







'6G1 


= —5u.Vu 


5u.Vq 




6G2 


= {K{qa) - K{q)) .Vqa 


5q.div Ua 


and < 


5G3 


= K{q).V6q 


q.div 5u 




6G4 


= {I{qa)- I{q))Aua 








= -I{q)A5u. 



Except Ra, all these additionnal terms are exactly the same as in [7], so we refer to this 
article for details on their estimates and we will only focus on what changes. As i] is 
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small, we can additionnally assume rj < 1, let us denote once again by C a constant only 

2 

depending on |^ and d (that may change from line to line). If we introduce for /i G [0, 1[ 

Ut)=^\\{6q,6u)\\ . 

eS (t) 

then using (2.10), we obtain that for all t G M (see [7] for details), 
fait) < Ce'^— " 

/ fa{T)F{T)dT+ [ ||<5n||.. (||g||.. , + (l + z.)||g||.d )(ir + ||i?„|| 
Jo Jo Baa \ ^2,1 ^2,1 ^t^L 



< Ce -0 



f Ur)F{T)dT + 7^-{l + -)Ut) + \\R^\\ 

Jo ^0 V L\Bl^ 



(3.81) 



where ^ 

F{t) = \\q\\ . . + ||n|| .d+i + (1 + -)\\ua\\ .d+i- 

^2,1 ^2,1 ^2.1 ^ ^2,1 

Remark 16 We recall that, due to endpoints in the Littlewood-Paley remainder term 
estimates, as in [7] we have the condition h < d — 1, which is why we impose h < 1 in 
order to work for any dimension d > 2. 

If 7] is so small that rjCe max(l, — )(1 + ^) ^ 2' ^^en we obtain thanks to the 

Gronwall lemma, 

Ut) < Ce^'^^'P.II . . ,_,e- "^^/o^W^^ 

Thanks to (3.78) and (3.79) we have, 

ft 



f F{T)dT< (- + -{2+-)] Ct] 
Jo \^ I'O ^ J 



so that using the condition on rj, we end up with: 

Ut)<C{^,^,d)\\R^\\ (3.82) 

Remark 17 For small viscosities, in the case A = 0, the previous condition on rj is 
roughly C-^e ^ < 2 which is obviously implied by the condition required in the existence 
result (see remark 15). 

To estimate the remainder in the case h g]0, 1[, we simply write that: 

which allows us to write: 

Hi? II H - K V2J'(f-'^-l)^^ll«-llr, - K V2J'(i+2)||r7.|U2^^^^ 

We now consider two cases: 
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If 2^ > a then 

If 2^' < a then as /i G [0, 1[ 

2i(2-^) 2J'(2-«) a^^-^") 

^ , = ^ < ^ < a"''. 

+ 2^-? 



We conclude that: 



^21 -62^1 



and then 

To complete the proof of the theorem, all that remains is to estimate 

Ca - p = * yOa - /9 = * (pa - p) + * /9 - /9 = (/"a * (^a - g) + * g - Q 

Thanks to what precedes only the last term has to be estimated, and similarly we have 



and 



D 2 o 2 

_Dr, , tin , 



For the hybrid norms, we do the same : 

Ua *q-q\\. . < Yl 2^^^-^-i) min(a2, 22^-) min(l, —)\\qj ^2 

= ^2^'(f+i)||g,-||i2a22-J'('^+2)min(l,^)2 < . (3.83) 

so that we finally obtain that if h g]0, 1[: 

Cyisc 

Ft? 4«; 

Coming back to (3.82) in the case /i = 0, we simply write that: 

22i 

and thanks to the Lebesgue theorem for series we obtain that the norm ||(/?q — p^Ca — 
p,Ua—u)\\ goes to zero when a goes to infinity, which ends the proof of the theorem. 
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4 Appendix 



The first part is devoted to a quick presentation of the Littlewood-Paley theory and 
specific properties for hybrid Besov norms used in this paper. The second section to 
general considerations on flows. 

4.1 Besov spaces 

4.1.1 Littlewood-Paley theory 

As usual, the Fourier transform of u with respect to the space variable will be denoted by 
J-{u) or u. In this section we will briefly state (as in [6]) classical definitions and properties 
concerning the homogeneous dyadic decomposition with respect to the Fourier variable. 
We will recall some classical results and we refer to [2] (Chapter 2) for proofs (and more 
general properties). 

To build the Littlewood-Paley decomposition, we need to fix a smooth radial function 
X supported in (for example) the ball i?(0, |), equal to 1 in a neighborhood of i?(0, |) and 
such that r i— )• x{''~-^r) is nonincreasing over M+. So that if we define ip{S,) = x(C/2) — x(C)' 
then if is compactly supported in the annulus G ^^^,00 = | < |^| < Cq = |} and we 
have that, 

VeeM'^\{0}, J^(/j(2"'0 = 1- (4.84) 

i&z 

Then we can define the dyadic blocks {Ai)i(zz by A; := ip{2~'-D) (that is A^n = 99(2~'^)2(,^)) 
so that, formally, we have 

?x = ^A,u (4.85) 

I 

As (4.84) is satisfied for ^ 7^ 0, the previous formal equality holds true for tempered 
distributions modulo polynomials. A way to avoid working modulo polynomials is to 
consider the set S'^ of tempered distributions u such that 

lim llSiullioo = 0, 

l—^—oo 

where Si stands for the low frequency cut-off defined by 5; := x(2^'-D). If u G 5^, (4.85) 
is true and we can write that Siu = A^n. We can now define the homogeneous 

k<l-l 

Besov spaces used in this article: 

Definition 3 For s G M and 1 < p,r < 00, we set 

if r < 00 and := sup2'*'||A;u||ip. 

p,oo ^ 

We then define the space B^ .^ as the subset of distributions u £ S'^ such that \\ 
finite. 
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Once more, we refer to [2] (chapter 2) for properties of the inhomogeneous and homoge- 
neous Besov spaces. 

In this paper, we mainly work with functions or distributions depending on both the 
time variable t and the space variable x. We denote by C{I;X) the set of continuous 
functions on / with values in X. For p E [1, oo], the notation X) stands for the set 
of measurable functions on I with values in X such that 1 1->- ||/(i)||x belongs to LP{I). 

In the case where / = [0,T], the space LP([0,T];X) (resp. C{[0,T]; X)) will also be 
denoted by Lj,X (resp. CtX). Finally, if / = M+ we alternately use the notation L^X. 

The Littlewood-Paley decomposition enables us to work with spectrally localized 
(hence smooth) functions rather than with rough objects. We naturally obtain bounds for 
each dyadic block in spaces of type L^L^. Going from those type of bounds to estimates 
in Lj.Bp^^ requires to perform a summation in ^^(Z). When doing so however, we do 
not bound the Lj,Bp.^ norm for the time integration has been performed before the £^ 
summation. This leads to the following notation (after J.-Y. Chemin and N. Lerner in 
[10]): 

Definition 4 For T > 0, s G M and 1 < r, p < oo, we set 

•= llr(Z)- 

It is then possible to define the space Lf^Bp ^ as the set of tempered distributions u over 
(0,r) X such that lim„^_oo SqU = in L''([0, T]; L°^(R'^)) and ||n||Tp < oo. The 

letter T is omitted for functions defined over The spaces Lj,Bp ,^ may be compared 
with the spaces Lj,Bp j. through the Minkowski inequality: we have 

||i*|l7PR3 <||'u||rPRs if r > p and ||n||rp >\\u\\rPn, if r < p. 

All the properties of continuity for the product and composition which are true in Besov 
spaces remain true in the above spaces. The time exponent just behaves according to 
Holder's inequality. 

Let us now recall a few nonlinear estimates in Besov spaces. Formally, any product 
of two distributions u and v may be decomposed into 

uv = TuV + TyU + R{u,v), where (4.86) 

TuV ■.= '^^Si-iuKiv, TyU :=''^^ Si^ivAiu and R{u,v) ■.= '^^ AiuAi'V. 

I I I |/'"«|<i 

The above operator T is called a "paraproduct" whereas R is called a "remainder". 
The decomposition (4.86) has been introduced by J.-M. Bony in [3]. We refer to [2] for 
properties, and also to [7] or [9] for paraproduct and remainder estimates for external 
force terms. 

4.1.2 Complements for hybrid Besov spaces 

As explained, in the study of the compressible Navier-Stokes system with data in critical 
spaces, the density fiuctuation has two distinct behaviours in low and high frequencies. 



35 



separated by a frequency threshold. This leads to the notion of hybrid Besov spaces and 
we refer to R. Danchin in [12] or [2] for general hybrid spaces . In this paper we only will 
use the following hybrid norms: 

=^ 5]min(a2,22i)2^-^||A,/||^2 
jez 

= ^ 2^-(«+2)||A,-/||i.+ «'2^-^||A,/||^.. (4.87) 

In this formulation, we obviously remark the threshold frequency log2 a which separates 
low (parabolically regularized) and high (damped with coefficient a^) frequencies. But 
as we prove in (2.30), the frequency transition is in fact continuous and the following 
equivalent formulations show that these norms are completely tailored to our capillary 
term: 

Exactly as in [9], the non-local capillary term a'^{(j)a*^Qa — ^Qa) has the same regularity 
as the capillary term from the local model V Ag: both of them belong to 1 1 ) • 



4.2 Estimates for the flow of a smooth vector-fleld 

In this section, we recall classical estimates for the flow of a smooth vector-field with 
bounded spatial derivatives. We refer to [15] or [6] for more details. We also refer to [22] 
for the incompressible Navier-Stokes case. 

Proposition 8 Let v be a smooth globally Lipschitz time dependent vector-field. Let 
Wit) := Jq \\S/v{t')\\Loo dt'. Let ijjt satisfy 

V't(x) = x+ [ v{t',4)t'{x))dt'. 
Jo 

Then for all t G M, the flow if^t is a smooth diffeomorphism over and one has ift > 0, 

\mt'\\L^<e^'^^ 

llZ^V^^i-Jdllioo <e^W-l, 

||I?V?1l«> < e^W ['\\D\{t')\\L^e'^^''Ut', 
Jo 

WD^^pt'h^ < e^(*) \\D\{t')U^e''^^''^dt' + 3(e^W \\DMt')\\L^e'^^''Ut' 

As in [6] we also introduce the two-parameter flow (t, t' ,x) 1— )■ X{t, t' ,x) which is (uniquely) 
defined by 



X{t,t',x)=x+ [ v{t",X{t",t',x))dt". (4.^ 
Jt' 
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Uniqueness for Ordinary Differential Equations entails that 

X{t, t", X{t",t', x)) = X{t, t',x). 
Hence ipt = X{t,0, •) and ip^^ = X{0,t,-)- 

Proposition 9 Under the previous notations, the jacohian determinant of X satisfies: 

det{DX{t,t',x)) = e//'(div^)(r,X{r/,x))dr^ (^gg^ 

and 

{ det{Dll}t{x)) = e/o(div«){r,i/..(x))dr^ 

\det{D'4)'[^{x)) = e" /o('^'^^)('^'"^(^'*'^))'^^ = e~ ^0^^"^^^'^''^'^°^*^''^^^'^'^ . 
Proof: differentiating (4.88) with respect to x, one gets by virtue of the chain rule, 

DX{t,t',x) = Id+ [ Dv{T,X{T,t',x)) ■ DX{T,t',x)dT. (4.90) 
Jt' 

This immediately implies that: 

dt{DX){t,t',x) = Dv{t,X{t,t',x)) ■ DX{t,t',x), 

and 

dtdet{DX{t,t',x)) = tr {Dv(t, X{t,t' ,x))) ■ det{DX{t,t\x)), 
so that we obtain the result. ■. 

4.3 Bessel functions 

In (2.7) we obtained that the interaction potential can be written the following way: 

'^(^) = , ^i, Kd^^{\x\), 

where K is the modified Bessel function of the second kind. In this section we will give 
specific properties for this Bessel function. 

Remark 18 We refer the reader to [1],[4], [27] or [36] (among a very rich litter ature) 
for a profusion of results and refinements for the Bessel functions. In this paper we will 
restrict to a very limited number of properties of function that will he used in section 
2.4.4. 

Let us begin with general facts on Bessel functions. If n is an integer, the Bessel function 
Jn represents the ?i-th Fourier coefficient of the function 9 i— )• e*^**'"^: 

Vx G M, Jn{x) = — [ e-^"^e^^^'°^d0 = - / cos{n9 -xsme)de. 
2vr Jo vr Jq 

ively, Jnix) can be sef 
X I— >• e^^*~"t^ into powers of t. 



Alternatively, Jn{x) can be seen as the coefficient of in the development of the function 
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Remark 19 This function can be extended into Ju{z) with z G C and for Re{u) > — ^ 
but in this paper we will restrict to real variable and an index G Z or | + Z. 

For a general index ly, J^, solves the following differential equation: 

x'^y"{x) + xy'{x) + (x^ - v'^)y{x) = 0. 

If ^ Z, a basis of the space of solutions of this differential equation is given by {J^, J-u)- 
If 1/ = n is an integer, we have J-n{x) = (— so that one introduced, as a second 
element of a basis of solutions, the following function (also called the Bessel function of 
the second kind) : 

^ _^ J.W -(-!)" 

u^n u — n 

In physics, many functions arise which are similarly constructed from the general Bessel 
function. Let us now consider the case of the modified Bessel functions of the first and 
second kind: for a general real index u > and a complex variable z we define 

/.(z)=e-^-V.(iz), and = ^^Z^lM^^lM. 

2, sm(z^7rj 

We can prove that (/^, Ky) is a basis of the space of solutions of the following differential 
equation: 

x'^y"{x) + xy'{x) - (x^ + v'^)y{x) = 0. 

We will now focus on Ky and recall some important properties (we refer to the cited 
books): 

Proposition 10 ([1],[4] or [36]) The function satisfies the following properties: 

2. For all x > 0, Ky{x) > 0. 

( — In X if u = 0, 

r(i^) /2\ 



4. K,{x) ~ . r/—. 

5. For all x ^W^, 

K,{x) = K,+2{x) - ^^^^i^.+i(^) = K,^2{x) - ^^^^K._i(x). 

X X 

6. Ky is decreasing from U*^ to M^: for all x £ M*^, Kq{x) = —Ki{x) and 

Klix) = -hKy_,{x) + Ky+i{x)) = -Ky^iix) - -Ky{x) = -Ky{x) - i^,+i(x). 
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7. For all X > and > 0, 



K,{x)= / 



-xcht 



ch{vt)dt, 



which implies that K^{x) is increasing with respect to v > (see [28]). 

We are now able to state the following estimates that will be needed in section 2.4.4. 
Proposition 11 Assume that d >2. 

1. There exists a constant Cd > such that for all x > we have: 



Kd_-^{x) > < 



Cd- 



X2' 



C- 



l 1 + 



ifd > 2, 



ifd = 2. 



2. For all r] g]0, 1[, there exists a constant Cd^r] > such that for all x > we have: 

g-(l-r;)a: ^-{l-r])x 

Kd_^{x) < Cd,n dT^, and \K'a__^{x)\ < C^_^ 

2 rjf, 2 2 J. 2 



Proof: For the first point, in the case = ^ — 1 > (d > 3) let us study the following 
positive function / : a; G i— t- f{x) = e^x^ Ky{x). Thanks to point 3 from proposition 
10, we have 



f{x) ^ r(z/)2'^-i > 0, 



and for the limit at infinity. 



( T 1 

OO \l U > ly, 



fix) ~ \ T^X" -> < 
x—^oo \ Z x—^oo 



ifi/<i. 



Using points 1,5,6 and 7 from proposition 10, we compute the derivative of /: which is 
positive : 



Vx > 0, f{x) = e'^x'^iK^x) - K^^i{x)) < 



'> \iv > \, 
= \iy = \, 
<0 \iv<\. 



Thus if 1/ > i, for all x > 0, f{x) > T{u)2''-^, that is when d>3, 



i^l_i(x)>r(^-i)2f-^4^, 



X2 



39 



Obviously when u = (that is d = 2) the previous study is useless and we need to 
consider g : x £ i— )• g{x) = e^(l + ^/x)Ko(x). For all x > 0, g{x) > and from the 
previous results, we have: 

g{x) +00, 
g{x) 



2 • 

So that either the lower bound m oi g satisfies m > either we have m < and 
then m is reached in some segment so that as g{x) > 0, we have m > 0. In any case we 
are sure that m > 0, which ends the proof of point 1. 

The proof of point 2 is also very elementary: for a fixed rj s]0, 1[, we introduce the 
function h : x £ M*^_ i— >• f{x) = e^^~'^^^x'^^'^Ky{x). For any > 0, we have: 

'h{x) ^0, 
h{x) 0, 

and h is then a bounded function, which implies the first estimate for u = f — 1- For the 
last point, writing: 

Kix) = ^K,{x)-K,+iix). 
and using what precedes, immediately implies that for all x > 0, 



and the proof is complete. 



The author wishes to thank Raphael Danchin, Nicolas Fournier, and Frangois Vi- 
gneron for useful discussions. 
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